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ABSTRACT 

We discuss the mass-radius (M-R) relations for low-mass (M < O.IMq) white dwarfs (WDs) of arbi- 
trary degeneracy and evolved (He, C, O) composition. We do so with both a simple analytical model 
and models calculated by integration of hydrostatic balance using a modern equation of state valid 
for fully ionized plasmas. The M-R plane is divided into three regions where either Coulomb physics, 
degenerate electrons or a classical gas dominate the WD structure. For a given M and central tem- 
perature, Tc, the M-R relation has two branches differentiated by the model's entropy content. We 
present the M-R relations for a sequence of constant entropy WDs of arbitrary degeneracy parameter- 
ized by M and Tc for pure He, C, and O. We discuss the applications of these models to the recently 
discovered accreting millisecond pulsars. We show the relationship between the orbital inclination 
for these binaries and the donor's composition and Tc- In particular we find from orbital inclination 
constraints that the probability XTE J1807-294 can accommodate a He donor is approximately 15% 
while for XTE J0929-304, it is approximately 35%. We argue that if the donors in ultracompact 
systems evolve adiabatically, there should be 60-160 more systems at orbital periods of 40 min than at 
orbital periods of 10 min, depending on the donor's composition. Tracks of our mass-radius relations 
for He, C, and O objects are available through the electronic version of this paper. 
Subject headings: binaries: close— pulsars: individual (XTE J0929-314, XTE J1751-305, XTE J1807- 
294) — white dwarfs — X-rays: binaries 



1. INTRODUCTION 

The discovery of three X-ray transient ultracom- 
pact accreting millisecond pulsars (MSPs), XTE 
J1751 -305 dM arkward TeTal] [200^ . XTE J0929- 
314 ( Gallowav ct al. 2002), and XTE J1807-294 
ijMarkwardt et al. 2003a b) have demonstrated the 
existence of binary pulsar systems with low mass, 
M2 w 10~^Mq, donors. These three ultracompact 
systems (here defined as binaries with orbital periods, 
Porh < 60 min) are remarkably homogeneous, with 
measured Porb — 42.4, 43.6, 40.1 min respectively, well 
below the minimum period for a syst em with a donor 
comp osed primarily of hydrogen (jRappaport et all 
Il982j) . Since the nature of the donors in these systems 
today depends on the prior evolution of the system, it is 
useful to discuss the potential formation mechanisms for 
these systems. 

Binary systems with Porb < 80 min can form 
through at least two channels. Stable mass trans- 
fer from an evolved main-sequence star (Nelson ^1;; ^1. 

t86; Fedorova & Ergma 1989; Podsiadlowsk i et al. 
02; Nelson & Rappa pord l2003() or a He burning star 
avoniie et al.. .1986.) onto a neutron star (NS) is one 
mechanism. In this scenario, the main-sequence star is 
brought into Roche lobe contact due to orbital angu- 
lar momentum losses from magnetic braking at a time 
when the core has nearly completed H burning. Such 
a system will evolve to orbital periods comparable to 
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the ultracompact MSPs and can reach Pnrh ~ 10 mi n 
jPodsiadlowski et al."2002'; 'Nelson & Rappaport"200%^ 
PodsiadlowskLct_al. (2002) and Nelson & Rappapo^ 
(2003) show that the resulting ultracompact binaries 
have donor masses M2 « 0.1 — O.2M0 as they pass 
through Porb ~ 40 min on their way towards a shorter 
period. These masses are significantly greater than those 
mea sured in the ult racompact MSPs ( Gallowav et ^ 
i002; Mar kwardt et al. 2002; Bildsten 2002). However, 
systems evolving through 40 min on the way out from the 
period minimum have masses more in line with the mea- 
surements (M2 ~ O.OIM0) and by this time the donors 
have become parti ally degenerate with core t emperatures 
Tc ~ 10^ - 10^ K ((Nelson fc Rappapo rt 2003). 

The second scenario that may form ultracompact 
systems involves triggering a common envelope phase 
during an unstable mass transfer episode from the 
donor onto the NS. The core of the donor, either a He 
or C/0 white dwarf (WD), and the NS spiral- in to 
shorter orbital p eriods until the envelope is expelled 
l)Paczvnskil I1976D . Several authors have proposed 
binary evolution scenarios in which the system, after 
emerging from the common envelope phase, then suffers 
in-spiral due to gravitational wa ve (GW) emission and 
eventually r e -estab l ishes contac t lllben fc Tutukovlll985t 
Rasio et al.l l2000t iDewi et all 120021 lYuneelson et alJ 
2002() . During this long GW in-spiral, the WD will 
have had time to cool, setting the entropy of the 
donor at the on set of the second mass transfer phase 
(|Bildstenl 12002(1 . iTaurisI ((1996) finds that a large 
fraction of the NS-WD binaries that undergo a com- 
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mon envelope phase will reach contact within 1 Gyr. 
Even considering a longer 4 Gyr delay between the 
formation of the WD secondary and the onset of mass 
tran sfer, a He WD will have T ,, k. 3 x 10^ - 10^ 
K lA lthaus fc Bcnvenuto "199/; "Driebe ct al.' Iggl 
iSerenelli. Althaus. Rohrmann, fc Benvenuto 2QQ1,1 . 
while a C/0 WD will have « 2 - 3 x 10^ K 
(jSalaris et alJl2000|) . The mass transfer time-scale at 
contact is much shorter than the WD cooling time-scale 
for these WDs so that the initial entropy of these objects 
is the minimum attainable. As noted bv iBildstc n (,2002), 
if these objects adiabatically expand under mass loss, 
their Tc will have been reduced by a factor of « 15 by 
the time they have reached Mi sa O.OIM0. 

In addition to the evolutionary arguments that donors 
in the ultracompact MSPs have not reached a T = 
configuration, the system XTE J1751-305 provides ob- 
servational evidence for a hot donor since, as noted by 
[Bildsten (2002), a fully degenerate companion composed 
of He or C can not fill its Roche Lobe (RL) in this sys- 
tem. Hence, in examining the donors in the ultracom- 
pact accreting MSPs, we need to consider them to be 
arbitrarily degenerate, low mass objects of evolved (He 
or C/0) composition. To further constrain the nature of 
these donors (for example, their Tc) requires knowledge 
of their mass-radius, AI-R, relation. For the composi- 
tions (He, C/0), mass (- O.OIMq) and (- 10^ - 10^ 
K) ranges of relevance to these objects, the correspond- 
ing central densities, pc lO'^gcm"^) are such that 
Coulomb and thermal contributions to the equation of 
state (EOS) provide non-negligible corrections to the de- 
generate electron pressure, impacting their M-R rela- 
tions. In this paper, we detail the M-R relation for low- 
mass stel lar objects of finite T r , exte nding the M-R re- 
lations of iZapolskv fc SalpeteJ l)1969D for T = objects. 
In particular, we make clear that there is a continuous 
connection between fully degenerate objects (i.e. WDs), 
fully convective low-mass stars (i.e. n = 3/2 polytropes), 
and Coulomb dominated objects. 

We begin in i)2.1l by constructing a simple model 
of these objects using an approximate EOS. Although 
crude, this model describes adequately the relevant 
physics and yields an analytic description of the qual- 
itative behavior of the M-R relations and how they 
are affected by Coulomb and thermal contributions 
to the EOS. We find that at finite Tc, arbitrar- 
ily degenerate sequences exhibit a two branch solu - 
tion, a fact noted pre v iously [e.g.lCox & Guili ( IQBS ) ; 
Cox fc Sabeteil (I1964D: [Hansen k Soangenber j ljl97Hi : 
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Tc, the sequence of solutions on these two branches ex- 
hibit a mutual end point at a non-zero mass, Mmin, 
below which equilibrium solutions do not exist. When 
Coulomb contributions are small and electrons are 
non-relativistic, fully convective stellar models of ar- 
bitrary deg eneracy are well represen t ed by n = 3/2 
polytropes ( Havashi & Nakano 1963"; 'Stevenson' 119911 

IBurrows fc Liebert 1993: Ushomirskv et al 199^ . In 

ii2.2[ we review the role played by degeneracy in 
determining the M-R relation for n = 3/2 poly- 
tropes and the existence of a two branch solu- 
tion for the polytrope M-R relation. Other au- 
thors have noted that for a given M, there is a 
maximum Tc that such polytrope models may have 



(iCox fc Guilil I196SI: iRappaport fc Josd 119841: iSteyensonI 
n99lUBurrows fc Lieberdll 99,'^ ITJshomirskv et alJII 99811. 

We connect the existence of a maximum Tc with that of 
Minin- In ^2.31 we construct realistic M-R relations us- 
ing an EOS for fully ionized plasmas. There we exhibit 
explicitly the impact of Coulomb physics on the M-R 
relations. Like the simplified and polytrope models, we 
find that the M-R relations of this model exhibit a two 
branch solution and a non-zero Mmin at high Tc. 

In §?? we apply our stellar models to the ultracom- 
pact MSP systems. For each of these systems, there is 
a donor (of some composition and Tc) which will fill the 
RL at the required Porb for any given orbital inclination. 
And, for a given composition, there is, in each system, 
a relation between orbital inclination and Tc. We exam- 
ine what constraints this places on the composition and 
Tc of the donors in these systems. For example, in XTE 
J1807-294, C and He donors can have any Tc, while an O 
donor's Tc has a minimum value. In XTE J1751-305, all 
He or C/0 donors must be hot. We also examine how the 
future evolution of these systems depends on donor com- 
position and Tc and highlight the fact that a multiple- 
valued M-R relation leads to a multiple-valued relation 
between Torb and the system's mass transfer rate, M . 
Finally, in the context of adiabatic evolution, we dis- 
cuss the expected number distribution of ultracompact 
systems as a function of Porb- In a steady state, this 
distribution depends almost solely on the response of 
the donor radius, P2; to mass loss through the quan- 
tity ripi = dlnP2/'^lnM2. The increased importance of 
Coulomb physics in C/0 donors alters tir as compared 
to He donors and the expected distribution for the two 
donor types differs dramatically. Depending on donor 
type, the relative number of systems at Porb ~ 40 min 
as compared to those at Porb ~ 10 min is « 60 for He 
donors and k, 160 for C/0 donors. We conclude in §?? 
by discussing future applications of and refinements to 
our models. 

2. MASS-RADIUS RELATIONS FOR LOW-MASS 
ARBITRARILY DEGENERATE STARS 

In the mass range of interest (M2 < O.IM©), M-R 
relations for various objects have previously been con- 
structed. For H rich objects. [Burrows et al.l ([20011) sum- 
marize the work that has been done on the structure 
and evolution of brown dwarfs and related objects. For 
these objects, the models include detailed treatment of 
the equation of state (EOS) and atmospheric physics. 
On the other hand, for objects with more evolved com- 
position (i.e. He or C/0) the theory is not so mature. 
Zapolskv & Saloeter (1969) calculated the M-R rela- 
tions for T = objects using the EO S they der ived from 
the T homas-Fermi-Dirac equation (.Saloeter fc Zapolskvl 
[19671) . This EOS treats in an approximate manner the 
corrections due to Coulomb interactions and exchange 
effects in a fully-degenerate plasma. Additional M-R re- 
lations pr oduced by several different T = EOS are pre- 
sented in lLai. Abrahams, fc Shap irol (Il99 1[). but do not 
differ appreciablv from the iSalpeter fc Zapolskvl ([19671 ) 
results. For the partially degenerate case, there is a large 
body of literatur e for He and C/0 core WDs more mas - 
sive than O.IMq fFontai ne et al.l2001l:[Panei et al. 20001. 
But, only recently have M-R relations for arbitrarily de- 
generate He and C/0 WDs with masses < O.IM0 been 
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calculated l)Bildstenll2002(l and only for a limited number 
of cases. Here we fill this gap by constructing low-mass 
WD models utilizing a realistic EOS. 

2.1. A Simple Model for Arbitrarily Degenerate Stars 

Degenerate stellar objects with 10"^ < M/Mq < 10"^ 
have central densities, pc ^ lO^-lO'^gcm"'^. In this 
density range, the relative energy contributions to the 
plasma from an ideal Fermi gas and Coulomb interac- 
tions can be comparable. Also, for these densities, at 
a temperature T ^ 10® K, the thermal contributions to 
the pressure are about 10% of those of the electrons. 
At lower densities. Coulomb and thermal effects become 
even more significant in calculating the EOS. To examine 
the interplay between these contributions and how they 
impact the structure of low-mass stars, we start with a 
simple EOS for a plasma composed of ions with charge 
number Z and atomic mass number A, 



P{p,T) = P,{p) + P,^{p,T) + Pcip) , 



(1) 



where Pe is the pressure of a fully degenerate non- 
interacting Fermi gas of electrons Pe — KeP^^"^, Pid 
is the pressure of an ideal gas of ions and electrons, 
Pid = KidpT, and Pc is the negative pressure contri- 
bution due to Coulomb interactions in the Wigner-Seitz 
approximation, Pc = —Kcp"^^^- Here = 3.323 x 
10^^(2//ie)^/'^ dyne cm'' g^^/'^ (/ie is the mean molecular 
weight per electron and equals A/Z in a single composi- 
tion plasma), Kid = 8.25xlO^(l+Z)/A dynecmg^^K^^ 
and Kc = 2.23 x IQ^"^ Z'^/^{2l peY'^Ay^ecm^ g"^'^ . 

Consider a one-zone stellar model, i.e. a spherical sys- 
tem characterized by a single pressure and density, P 
and p. From dimensional analysis, P ~ GM'^/R^ and 
p ~ M/R^ where M and it! are the mass and radius of 
the star and G is the Newtonian constant of gravitation. 
With the pressure given by equation 



G^^K,dT—+K,-r^-K 



c- 



(2) 



i?5 i?4 

In the r limit, the term involving i^id vanishes and 
solving for R gives 



R{T = 0) = 



(3) 



Kc + GAf2/3 

showing that the gravitational and Coulomb attraction 
act to collapse the star, which is supported by the de- 
generate electron pressure. As M ^ 0, Coulomb forces 
dominate gravity and R cx A/^/'^. As M — + cx3. Coulomb 
forces become negligible and R cx The transition 

from gravitational to Coulomb dominated regime occurs 
where M ~ [Kc/Gf/"^ - 1Q-'^Z{2/ pefMQ. The rela- 
tion R oc M^/^ in the low M limit implies a constant 
density, Pmin, fixed by the balance between the electron 
pressure and Coulomb attraction. Exhibiting p explic- 
itly as a function of M, p = (GA/^/^ -1- Kc)lKe and 
Pmin = [Kc/Kef = 0.671Z2 gcm-3. 
For T > 0, equation (jSJ has solutions given by 

R^^i!!l(Kc + GM'^^± 

(4) 



2KidT 



AKeKidT) + 2KcGA'P/^ 



GH'I'^/^ ) 



o 
K 

o -1 




-6 -4 -2 

log(M/Mj 

Fig. 1. — The schematic M-R relation of our simple model given 
by equation jl}. The curves show isotherms at T = 10'', 10^, and 
10*' K for pure He (solid) and C (dashed) compositions. The curves 
on the upper branch are labeled with their respective temperatures. 
The lower branch curves at low M have a T = 10^ K. 



Figure n exhibits several isotherms of the AI-R relation 
of equation (01 . The two-branch nature of this relation 
is apparent, as is the clear separation into three classes 
of stellar objects. The large M lower-branch is made up 
of gravitationally bound objects supported by degener- 
acy pressure. On the lower-branch at small M, we have 
Coulomb dominated objects. The upper branch consists 
of objects supported by thermal pressure. 

For r > 0, there exists a minimum mass, Afmin, found 
by setting the discriminant in equation Q to zero. 



2{KeK,dTf'^ - Kc 
G 



(5) 



This expression is only positive if T > 4.53 x 
lQ^AZ^/^/{\ + Z){2/pe) K. For He(C) this temperature 
is 1.5 x 10^ (8.5 x 10^) K. Above this critical temperature, 
the two branches meet at a mutual end point and no so- 
lution with M < Afmin exists. When Coulomb physics 
is negligible, A^min occurs at the point where Pe — Pid 
along the solution curve. The existence of Afmin results 
from the fact that for Af < A/min, the pressure provided 
by Pe + Pid at any p is greater than that needed to sup- 
port the star (this is not the case for either an ideal gas or 
a Fermi gas independent of the other — in either of these 
cases it is well known that equilibrium solutions down to 
A/ = exist). Alternately, for a given A/, the isotherm 
on which M = Afmin gives the maximum T for which 
solutions exist with this mass. For an object starting 
out on an upper branch solution, e.g. a recently expired 
star, as it loses entropy via radiation, it contracts. For 
a fixed mass, a star supported by thermal pressure has 
T (X and it heats up, as expected. The thermal 

pressure in this case goes as Pid oc P~'*, but Pe oc R~^ 
and eventually Pe dominates, halting substantial con- 
traction. Further entropy loss leads to a reduction in 
T. For a given T and Af , the two possible solutions are 
physically distinguished by their entropy. 
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2.2. Polytrope Models for Low-Mass Stars Neglecting 
Coulomb Physics 

Ignoring non-ideal effects and assuming that the elec- 
trons are non-relativistic, a fully convective object has 
an EOS obeying P oc p^/"^ and is modeled by an 
n = 3/2 polytrope. The specific entropy, s, through- 
out such a model is a constant. As s only depends on 
the degeneracy parameter rj, defined as the ratio of the 
chemical potential of the electrons to kT (k being the 
Boltzmann constant), t? is also a con stant of the model 
llUshqmirskv et aTl Il998j) . Following lUshomirskv et aP 
l) 1998ft . this allows us to write the pressure of a non- 
interacting gas of electrons and ions as 

P 



P^—t—kT, (6) 

/ioffTOp 

where nip is the mass of a proton and /icft is defined as 



(7) 



J_ _ 1 ^ 2F3/2(??) 

thus varying as rj changes. Here /i^ is the ion mean 
molecular wei ght and Fk i s the F ermi-Dirac function of 
order k from iCox fc GiiilJ l)1968D . We then utihze the 
n = 3/2 polytrope relations, equation and the fact 

that Fi/2{ri) p/ ip-eT'^^'^) to determine M and i? as a 
function of r/, Tc- These results can be expressed as 

-1/3 



R 

Rq 



0.359 



M 



O.OIM, 



Tr = 3.42 X 10^ K 



M 



4/3 



2/3 



PcsF^ 



2/3 
/2 



O.OIMq^ 
rewrite 



1/3 



-,1/3 



the results 



(8) 



(9) 



which is a minor 
lUshomirskv et al.l l|1998 

Equation @ shows the relation between Tc and M is 

1/3 

a function of rj through the combination ^csF^^^ ; which 
has a single maximum at 77 w 3 — 5 for expected WD in- 
terior compositions. Just as in the simple models of t)2.1l 
this shows explicitly the connection between a maximum 
Tc for a given M and the existence of an Mmin for a fixed 
Tc- Here again, there are 2 equilibrium radii for a given 
M and Tc differentiated by their degree of degeneracy or 
equivalently their entropy. 

The transition from a thermal pressure dominated to 
degeneracy pressure dominated state in these models ap- 
proximately determines M„iin- This transition occurs 
near where « 2F3/2/3/Xei^i/2- For increasing fii, 

this occurs at lower rj, i.e. at lower pc if fJ^e and Tc are 

held fixed. With P cx p^^^, dimensional analysis shows 
1/2 

M oc Pc ; a lower density at the transition between de- 
generate and non-degenerate states gives a smaller Mmin. 
For a fixed Tc, a pure C WD will has a smaller M^i^ than 
a pure He WD. 

2.3. Mass-Radius Relations for Isentropic White 
Dwarfs 

We now calculate the M-R relations derived from stel- 
lar models that include Coulomb physics. In these cal- 
culations, we assume that the interior profile of our stel- 
lar models are adiabatic. In reality, the actual entropy 
profile in a given donor will depend on its evolution- 
ary history. Many factors — such as whether the system 



formed through a stable or unstable mass transfer chan- 
nel, whether or not H burning is still ongoing at the point 
of contact, and how the mass transfer rate changes with 
time — can impact either the initial entropy profile of the 
donor or its subsequent evolution. In calculating mod- 
els for donors in ultracompact binaries, without choosing 
their evolutionary history, a reasonable approximation of 
their internal entropy profile is the best that can be done. 
Since we aim to construct models that will enable anal- 
ysis of the donor's properties today, irrespective of their 
past histories, we must assume some internal profile in 
calculating them. 

A system that initiates mass transfer at a Porb ^ 40 
mill has a donor mass M2 > O.OIMq. The mass trans- 
fer time-scale for such a system with a NS primary is 
roughly M2/M < 1 - 100 Myr, depending on M2. Con- 
sideration of the flux through the half-mass point in our 
models due to electron conduction (calculated using the 
conductive opacities of'Potekhin et al. (1999)) compared 
to the heat content of the interior half of the model shows 
roughly that the time required to transport this heat out 
of the interior is ^ 100 Myr-1 Gyr, again depending on 
Al2. Thus during the mass transfer episode that would 
lead up to the creation of systems at Porb ~ 40 min, the 
internal evolution is to first order an adiabatic expansion 
and, in the absence of tidal heating, the initial entropy 
profile should be more or less preserved with some cor- 
rections due to heat transport. To what ever degree this 
occurs, the critical point is that the interior will not be 
able to maintain an isothermal profile. We chose to use 
an adiabatic profile, instead of another arbitrary choice, 
since it provides several convenient features — it is a lim- 
iting case for the possible thermal profiles of the donor 
and produces the most compact configuration for a given 
Tc and M2, it is completely determined by the utilized 
EOS, and it allows parameterization of a set of models 
by one quantity, the specific entropy which is constant 
throughout the model. In addition, the calculations of 
Nelson & Rappaport ( 2003) show that He donors tend 
to become adiabatic as they lose mass. These latter cal- 
culations also highlight that throughout the mass loss 
episode, donors remain far from isothermal except in the 
deep int erior of the star, and then only at masses near 
0.1 Mq ijNelson fc Rappaportil200^ . 

2.3.1. Equation of State 

For the equation of sta te, w e use the results of 

t.abrier & Potekhin (1998) and iPotekhin fc ChabrieJ 
)00). Their work provides accurate prescriptions for 
calculating the EOS for a fully ionized plasma in either 
liquid or solid states and which includes the ideal contri- 
butions from non-degenerate ions, degenerate electrons 
and the non-ideal contributions due to Coulomb inter- 
actions. We also include the radiation contribution, a 
small effect. In the calculations below, we assume that 
the plasma is in the solid state when the quantity 



akl 



(10) 



ijFarouki fc Hamaeuchil [T993f) . Here a is the inter- ion 
spacing given by a = (47mi/3)~^/'^, being the num- 
ber density of ions. The quantity F is a measure of the 
strength of the ionic Coulomb interactions. We display 
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where r2 is the adiabatic exponent 



d InP 
d\nT 



-1 1 

log(p/g cm-3) 

Fig. 2. — The He EOS we utilize showing the comparison 
isotherms (solid/dotted lines) and a set of representative adiabats 
(dashed lines). The isotherms are for temperatures incremented 
by Alog(T/K) = 0.5, with the upper curves corresponding to 
log{T/_R') = 7. Along the isotherms, the dotted lines indicate re- 
gions where the plasma is not fully ionized. The adiabats typically 
cross into regions where full ionization of the plasma is not definite 
and is a source of uncertainty in some of our models. 



a set of isotherms for a pure He plasma calculated with 
this EOS in Figure|21with the solid and dotted lines. The 
dotted lines indicate regions in parameter space where T 
and p are such that the plasma will likely become only 
partially ionized. In these regimes this EOS is not strictly 
valid. 

The dashed lines in Figure [5] show representative adi- 
abats for this EOS, some of which cross into regions of 
partial ionization. Since wc utilize an adiabatic inter- 
nal profile in calculating our models in t)2.3.2l in some 
of these models there can be a point in their outer lay- 
ers where our EOS becomes invalid. We use our EOS in 
these regimes despite this problem for two reasons. First, 
in the models of most interest only the very most outer 
layers of the models lie in regimes where ionization state 
transitions become an issue and the calculated A4-R rela- 
tion is not affected. Second, this EOS provides a simple 
method of calculating the EOS for various compositions 
over a wide range of p and T and this ease of use would 
be sacrificed by constructing composition specific exten- 
sions to the EOS. We will highlight in our results models 
in which these EOS concerns may cause more than a few 
percent uncertainty in the calculated M-R relations. 

2.3.2. Calculation of the Models and Results 

We constructed models for arbitrarily degenerate ob- 
jects by integrating mass conservation and hydrostatic 
balance while presuming an adiabatic temperature pro- 
file. The calculation of our models proceeds as follows. 
In the interior, where degeneracy pressure dominates, we 
integrate 

dT T2-lTdP ^^^^ 



(13) 



dr 



P dr 



found from the EOS. The density at each integration 
step is then solved for from P and T using the EOS. 
At low p, where Coulomb effects on the pressure become 
significant, we switch to integrating 

dp _ 1 p dP 

dr Ti P dr ' 

where 

is another adiabatic exponent, and solve for T at each 
integration step from P and p. Each model integration 
was terminated once the following criteria were met — 
(1) between two integration steps in which P differed by 
more than a factor of ~ 20%, m and r changed by no 
more than 1 part in 10^, and (2) the current pressure is 
such that P/Pc < 10"^ 

We change the integrated variable because as P ^ 
along an adiabat, p becomes very insensitive to P and it 
is numerically intractable to determine p by root finding 
from P and T. On the other hand, in the low P limit, 
determining T accurately from P and p is possible, some- 
thing that was not true in the highly degenerate regime. 
We switch from integrating T to integrating p when a 
rough measure of the degeneracy, 167 p/ p-eT^^^ < 100, is 
satisfied (where p and T are in cgs units). In the solid 
state, when F 3> 180, F2 ^ 1 due to the rapid decline 
in the plasma's specific heat once crystalline. We use 
an adaptive step-size explicit Runge-Kutta algorithm to 
integrate our equations with the step size chosen to main- 
tain a fixed fractional accuracy. A dramatic increase in 
F2 causes a sharp decline in the speed of the integration 
as the algorithm tries to maintain this accuracy in all 
three integrated quantities, to, P, and T. To deal with 
this problem, once T = 100 K, we set dT/dr « 0. This 
causes no issues in the M-R relations because by this 
point we are already well in the T — > limit as far as the 
P-p relation is concerned in any of our models. 

Finally in these models the ion coupling parameter, 
F, increases with r. This is due to the fact that along 
an adiabat, F3 — 1 = {d\nT/dlnp)ad varies from 2/3 to 
« 1/2 in a Coulomb plasma. Since F oc p^^^/T, and 
T oc p^^~^, F goes as p to a negative power. This is not 
just true for adiabatic profiles — any object with a profile 
T (X p^ for 7 > 1/3 will have F increasing with r (assum- 
ing Z is fixed) — because of this some models transition 
from liquid to solid in their outer layers. For these cases 
we do not attempt to match adiabats in the solid and 
liquid phase (i.e. we do not account for the latent heat). 
Instead, the integration assumes continuity of P and T; 
the entropy in these models have a small discontinuity. If 
crystallization of the object were to actually occur from 
the surface inward, this could have significant impact on 
the mass loss rate, since the primary's gravitational field 
would have to overcome the Coulomb binding of the crys- 
tal to effect mass transfer. More realistic calculations are 
obviously needed to consider this further and considera- 
tion of this potential effect in evolutionary calculations 
are encouraged. 
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Fig. 3. — A Comparison between M-R relations for polytropes 
and the full adiabatic models with Coulomb physics for pure He. 
The short dashed curves show the M-R relation for n = 3/2 poly- 
tropes of pure He at Tc = 10'', 10^, and 10^ K . The solid curves 
show He WDs calculated with the full EOS of i|2.3.1l at the same 
set of Tc as the polytropes. The significance of Coulomb interac- 
tions at the masses shown on the M-R relation is obvious. The 
dotted portions of the full curves indicate models where more than 
5% of the model's mass is located in regions where the EOS is not 
strictly valid. 




log(M/Mj 

Fig. 4. — The M-R relations for adiabatic models composed of 
He (solid lines) and C (dash lines). The dotted portions of these 
curves indicate models where more than 5% of the model by mass 
is is located in regions where the EOS is not strictly valid. The He 
curves show models with Tc = 10^, lO'^, 5 X lO'^, 10*^ and 5 X 10*^ 
K. The C curves, models with Tc = 10*, 10'^, 3 X 10^, and 5 X lO" 
K. The dot-dash lines are the M — R relations for objects filling 
their Roche lobes at the noted orbital period. 



Typical results for our model calculations are shown 
in Figures 13 and 0] Figure |31 shows for pure He models 
how the results utilizing the EOS of i )2.3.1l differ from the 
polytropes that neglect Coulomb physics. In this, and in 
Figure 01 portions of the M-R curves shown in dotted 
lines indicate models where more than 5% of the mass 
lies in regions where our EOS is not strictly valid. The 
impact of Coulomb interactions on the structure of low 
mass WDs is clear from the comparison of the He poly- 
trope models (short dashed lines) to the realistic EOS 
He models (solid lines). Figure 01 displays a set of M-R 
isotherms for our He (solid lines) and C (dash lines) mod- 
els along with lines of constant Porb for a donor filling its 
Roche lobe overlaid (dash-dot lines with Porb indicated) . 
Again, the dotted portions of the M-R relations indicate 
models in which more than 5% of the model's mass lies 
in a regime where the EOS is not strictly valid. From 
Figures and 01 it is clear that the realistic treatment 
of Coulomb physics in the EOS is necessary to calculate 
accurately the structure of low mass WDs. 

3. APPLICATION OF THE MODELS TO ULTRACOMPACT 
ACCRETING MILLISECOND PULSARS 

We now apply the adiabatic models of W2.3\ to the 
three known ultracompact accreting millisecond pulsars 
and the high-field X-ray pulsar 4U 1626-67. In Figure 
we display (short dashed lines) the AI-R relations of 
Roche lobe (RL) filling donors in XTE J0929-314 XTE 
J1751-305, and XTE J1807-294 j Markw ardt et alil200l 
iGallowav et~al]l2002l:lMarkwardt et al.ll2003albD and M- 
R relations for our adiabatic models. We show two 
isotherms each for He, C, and O models, an approxi- 
mate T = M-R relation and one for hot (T^ = 3 x 10^ 
K) models. If the donors in all three systems are He 



WDs, then T = objects are allowed in XTE J0929-314 
and XTE J1807-294, but XTE J1751-305 requires a hot 
donor (Bildsten 2003). For C/0 donors, whose M-R re- 
lations will lie between the C and O models shown, only 
XTE J1807-294 permits a RL fiUing cold donor. The 
other two systems both require hot C/0 donors. The 
curves showing the RL filling AI-R relations are param- 
eterized by the orbital inclination, i (where « = is a 
face on system) and there is a correspondence between 
Tc and i for each donor composition, as shown in Figure 
From Figure jH the donor in XTE J 175 1-305 must 
have Tc > 10® K. From orbital inclination constraints, 
the probability that XTE J1807-294 can accommodate a 
He donor is 15% (for XTE J0929-304, it is 35%). 

The other ultracompact ac creting pulsar, 4U 1626- 
67, at Pnrh =^ 41-4 min l|Middl editch et al.l 119811: 
lChakrabartvlll99^ is also shown in Figure |S1 Though 
the orbit has not yet been detected by timi ng the pulsar, 
the c u rrent upper limit of sin i < 8 It-ms l)Levine et alJ 
Il988t IChakrabartv et all Il997|) allows us, in conjunc- 
tion with our theoretical work, to constrain the nature 
of the donor s tar. Ever since the discovery of a neon 
emission line ()Angelini et al.l Il995() from this system, 
there have been active d i scussi ons of the nature of the 
donor. The iSchulz et al.l l)2001() measurement of a high 
Ne to O ratio (further inferred in other ultracompacts by 
Duett. Psaltis. & Chakrabartv (2001)), led them to sug- 
gest that the donors in these binar ies are the cores o f 
previously crystallized C/0 WDs. iHomer et al.l l|20fl2f) 
have since also seen strong C and O lines, but no evi- 
dence for helium. Hence, this system seems a likely one 
to use for probes of C/0 donors. 

For any sini, a star that fills the RL at the measured 
Porb for each of the ultracompacts can always be found by 
some combination of entropy, composition, and mass and 
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Fig. 5. — A comparison between the M-R relations for our 
adiabatic He (solid lines), C (long dashed lines), and O (dash-dot 
lines) along with M-R relations for Roche lobe filling donors in the 
three known accreting millisecond pulsars and 4U 1626-67 (short 
dashed lines). The He M-R relations have Tc = 10^ and 3 X 10^ K; 
The C and O, Tc = 10'* and 3 X lO*' K. The low Tc relations show 
approximately the T = relation for each composition. The RL 
filling solution for 4U 1626-67 extends down to Af = to indicate 
that this system has not yet had its mass function measured. 
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Fig. 6. — The relation between the orbital inclination, i, and Tc 
of our adiabatic He, C, and O donors in the ultracompact systems. 
Hot donors are required in XTE J1751-305 for either He or C/O 
WD donors; C/O donors in XTE J0929-314 must be hot, while 
for XTE J1807-294 a T = C/O donor is allowed if the system is 
nearly edge on. 



the current values of each impact the orbital evolution of 
the system. For M2/M1 < 0.8, the Roc he radius, 
can be approximated by l|Paczvnskilll967|) 



Rl w 0.46a 



1/3 



(15) 



where a is the separation between Mi and Af2- Combined 
with Kepler's third law, this leads to the so-called period- 



8.9 hr 



Rr. 



Mo 



(16) 



Assuming conservative mass transfer, the mass tran sfer 
rate (a positive quantity) is given by (IVerbuntlll99l) 



M 1 

M2 J riR, 5/3 - 2M2/M1 



(17) 



where J is the orbital angular momentum, J the 
angular momentum los s rate set by GW emission 
ijLandau fc Lifshit jll962() and 



d lni?2 
din Mo 



(18) 



denotes how the donor's radius changes with mass loss. 

For a given system, J will depend on the inclination 
through M2 and a. The rate at which the orbit evolves 
will vary accordingly as will the M—Porh relation over the 
course of the evolution. To illustrate this, we calculated 
the forward evolution of XTE J0929-304 assuming a He 
donor and four different sini values. We assumed that 
the NS has Mi = 1.4Mq (and ignore the change in Mi 
due to accretion) and set 



d lni?2 
dlnM2 



(19) 



ad 



SO that the donor evolves adiabatically, ignoring any 
heating mechanisms (such as irradiation or tidal heating) 
and cooling. We show the results in Figure displaying 
M2 and M as a function of Porb • These relations are not 
single valued, but parameterized by orbital inclination, 
or equivalently, by the donor's entropy. A smaller sini 
requires a more massive donor (which must be hotter 
than a lower mass donor if it is to have the fixed mean 
density imphed by the system's Porb; see also Figure EJ. 
This gives a higher M for a fixed Porb , as seen in Figure 
[7| This also impacts the rate at which the orbit evolves. 
In Figurc[3 the age of the system from today is indicated 
along each curve by symbols and it can be seen that the 
smaller sin i the faster the system will evolve in Porb ■ 

At a given orbital inclination in a specified system, a 
C/O donor must have a higher Tc to fills the Roche lobe 
than a He donor. This is due to the stronger Coulomb 
physics in the C/O object, which also causes riad to dif- 
fer between the donor types and impacts the binary's 
evolution. The difference in riad between composition is 
shown in Figure|S|for two representative adiabatic tracks 
for each composition. The solid dots show the Zapolsky- 
Salpeter for the same compositions. The effect of 
different riad on orbital evolution is evident in Figure 
El which compares the evolution of XTE J0929-304 for 
sin I = 0.6 and He, C, and O donors. The difference in 
the tracks comes about due to the difference in rzad due 
to composition; higher Z donors evolve fastest in mass, 
but slowest in Porb because they remain more dense than 
lower Z donors at a given mass. 

4. THE PERIOD DISTRIBUTION FOR ADIABATICALLY 
EVOLVING ULTRACOMPACT SYSTEMS 

These evolution calculations highlight the dependence 
of observables (Porb and M) on the donor's entropy 
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Fig. 7.— The relation between M2, M and Poih for XTE J0929- 
314 over the course of its orbital evolution forward from today as 
it depends on sini, or equivalently, the entropy of the donor. In 
particular, the M — Porb relation is parameterized by the donor's 
entropy. Along each curve, the symbols designate time from today: 
triangles, 500 Myr; squares, 1 Gyr; pentagons, 5 Gyr; circles, 10 
Gyr. A smaller sin i requires a more massive donor and will produce 
a faster evolution in Poib for the system. 
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Fig. 9. — A comparison between the evolution of XTE J0929- 
304 assuming different donor types for an orbital inclination of 
sini = 0.6. Shown are M2 and M as a function of Porb for He, 
C, and O donors. The contrast in the initial evolution between 
the three donor types comes from the dilferences in their n^d ■ The 
symbols mark the age of the system from today. 
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Fig. 8. — The adiabatic change in donor radius with respect to 
mass, Jiad for He (solid lines), C (dashed lines), and O (dashed- 
dotted lines) stars along two representative adiabats for each com- 
position. The lower curve in each set has Tc = lO'^ K at the high 
mass end; the upper curve has Tc = 10^ K (i.e. it is elfectivcly 
the T = sequence for each composition). The stronger Coulomb 
physics decreases the magnitude of the radius response to mass loss 
in the C/O WDs as compared to He WDs. The large dots show 
for the T = Zapolsky-Salpeter models for He, C, and O (bottom 
to top). 



and composition. Wc now emphasize the impact of the 
donor composition on the ultracompact population, es- 
pecially the resulting orbital period distribution. In the 
scen ario of cooling; WDs re aching contact via in-spiral 
(e.g. iNelemans et al.l l|2001|) '). the relative number of He 
vs C/O WD's that come into contact and stably reach 



longer orbital periods is hard to know. However, what we 
will show here is the ability to constrain the relative pop- 
ulations of say, He donors, at one orbital period versus 
another. 

A more complete picture of this dependence is shown 
in Figure [TUI where we display for He (solid lines), C 
(dashed lines), and O (dash-dot lines) donors, the M- 
Porb relation along M-R isotherms with Tc = 10^, 3 x 
10^, 10'^, and 3 x 10^ K, assuming Mi = 1.4 Mq and n = 
Had- These are instantaneous M along an adiabatic track 
at a given Porb and Tc- One can see immediately that for 
a given Porb, M can constrain the combination of donor 
Tc and composition. In particular, a sufficiently strong 
upper limit on M can rule out a He donor for a given 
Porb- Above the minimum M for a He donor, further 
information about the donor composition is difhcult to 
infer without constraints on Tc- 

Now consider adiabatic evolution with initial Mi — 
1.4 Mq and donors of varying composition and Tc which 
fill their RL at Porb = 10 min. We evolve these systems 
assuming the donor responds adiabatically to mass loss. 
The resulting tracks in the M-Porb diagram are shown 
in Figure 1111 along with the measured periods of the 
known ultracompact binaries with a NS primary (verti- 
cal dotted lines) and the critical M below which the ac- 
cretion disk in these systems is subject to t hermal insta- 
bilitie s for both irradiated (hatched region, IDubus et alJ 
^9991)) and non-irradiated disks (nearly horizontal lines, 
iMenou et al (2002)). As compared with Figure [TOI it 
can be seen that for systems at Porb > 30 min to have 
time-averaged mass transfer rates M > 1O"^°M0 yr'^ 
the donor cannot have evolved adiabatically from sys- 
tems coming into contact at Porb ~ 10 min. There 
are potentially two examples of such systems: 4U 1626- 
67 (Porb = 41.4 min and M > 2 x Mq yr''^ 
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Fig. 10.— The M 

"-forb relations assuming npt — nad for He (solid 
lines), C (dash lines), and O (dash-dot lines). For each composi- 
tion, the four curves show lines of constant Tc = 10^, 3 X 10®, 10^, 
and 3 X 10'^ K (bottom to top). The bolder lines are the T = 10^ 
K curves. 



llChakrabartv et al.llT997l) ) and 4U 1916-05 (Ppi-b = 49.8 
min, M 5 X 10''^° Mqyt-^ ijSwank. Taam. fc White! 
119841) ). If these measured Ms reflect the long-term av- 
erage M, then the donors must be extremely hot (> 10^ 
K, see Figure However, the typical M at these 

orbital periods are below where a He or C/0 disk is 
thermally st a ble (s ee Figure ^ iMenou etall l)2002D : 
iDubus et al.1 ((1993 )). In that case we would explain 
the present M as a higher rate indicative of a system 
in outburst. Indeed, the luminosity from 4U 1626-67 is 
observed to be in a steady slo w decline ijMavromatakisI 
IT99I IChakrabartv et alj[l997t) . 

If the ultracompacts evolve adiabatically, we can deter- 
mine their relative numbers as a function of Porb- Defin- 
ing iV(Porb) such that N dPorh is the number of systems 
with orbital period between Porb and Porb + rfPorb and 
demanding continuity gives 



d{NP, 



orb ) 



dR 



= 0, 



orb 



(20) 



leading to the expected relation between Porb and N: 



N 



Po 



b,0 



Pn, 



(211 



where A*o and PorbiO are the respective quantities at some 
reference orbital period, Porb,o- For nad fixed, R cx M"°<i 
and with M2 < Mi, equations and |(2U lead 

to the simple relation 



d\nN _dln Porb 
d In Porb 



11/3 - 5nad 



d In Porb 
In this case 



1 - 3n, 



ad 



N 



p 



orb 



(22) 



(23) 



,Porb,0, 

and from this it is clear that rtad alone determines the 
number distribution. As a increases with riad, systems 
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Fig. 11.— The M 

"-forb relations along adiabatic evolution- 
ary tracks. Each track starts with a donor filling its Roche lobe 
at Porb = to min. For each composition — He (solid lines), C 
(dash lines), and O (dash-dot lines) — tracks for models with initial 
Tc = 10^, 3 X 10^, 2 X lO**, 5 X 10® and lO'^ K are shown (bottom 
to top). By the time these systems have evolved to Porb > 30 
min. A/ < IO^^^A/q yr~^; donors in binaries with Porb ^ 30 min 
that have persistent M's higher than this cannot have adiabati- 
cally evolved from Porb < 10 min. The vertical dotted lines show 
the orbital periods for ultracompact systems with a NS primary. 
For each composition, the upward sloping lines show the critical M 
below which the accretion disk is thermally unstable ignoring irra- 
diation I Mcnou et al 2 00^ . The shaded horizontal band gives the 
critical M for an irradiated disk IDubus et al.lT999 ) for a range of 
irradiation efiiciencies. The ba nd corresponds to the range of val- 
ues for the IDubus et alj tT999l) C parameter of ±50% the fiducial 
value. 



with C/0 donors will have a stronger increase in TV with 
Porb than those with He donors due to the difference in 
riad shown in Figure |S1 In the more general case, nad is 
variable and equation H22() becomes, up to terms of order 
M1/M2, 



11/3 - 5nad 



18 



dnad 



dlnN 

din Porb 1 - 3nad ' (3nad + 5)(3nad - 1) In Porb ' 

(24) 

and in general the distribution is almost solely a function 

of nad- 

While equations H22|l and (|24|l highlight the central- 
ity of nad in determining N, to calculate N for each of 
the adiabats in Figure ITTl it is more straightforward, and 
slightly more accurate due to the small dependence of N 
on M2/M1, to calculate Porb numerically and evaluate 
equation (|21|l . We calculate Porb from 



Porb = — 



dR 



orb 



dAh 



Mo = 



P 

'2M2 



(3nad-l)M2 



(25) 



where M2 = — M. We display the resulting N distribu- 
tion for each adiabatic track, normalized to the value of 
A^o at PorbjO — 10 min, in Figure El It is clear that the 
distribution is a strong function of composition, through 
the differing n^ = nad, but only a weakly depends on the 
entropy of the donor. The entropy dependence of N de- 
rives from the fact that Porb oc M2 [equations H15|l . (|17|l . 
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Fig. 12. — The number distribution, N{Poih) of ultracompact 
systems along the adiabats shown in Figure [TTl Each distribution 
is normaUzed to 1 at Porb = 10 min. For each composition, from 
top to bottom, the tracks are for donors with initial Tc = lO'^, 3 X 
10'', 2 X 10^, 5 X 10^ and lO'' K. The differences in the distributions 
between donor types are caused by the differences in their riad- 
The differences in n^d can be seen to play a more significant role 
in determining N{Parh) than that of the initial donor entropy. 



and H25|) ] and donors with a higher entropy at a fixed 
mean density have larger M2 (see Figure ISJ. Figure IT^ 
also shows that depending on donor type and entropy, we 
expect to see roughly 60-160 (for He and C/0 donors, re- 
spectively) times as many systems at Porb ~ 40 min than 
at 10 min. While consideration of non-adiabatic evolu- 
tion will change the value of this ratio, the fact that it will 
be larger for C/0 donors than for He donors is expected 
to be a robust result, regardless of the evolutionary sce- 
nario considered. 

5. DISCUSSION AND CONCLUSIONS 

We have presented models for arbitrarily degenerate 
stellar objects including Coulomb physics with masses 
M < O.IM0. At these low masses, the well known 
AI-R relations for n — 3/2 polytropes (i? oc M), 
WDs (i? oc Af^^/'"^), and Coulomb dominated objects 
(i? oc M^/'^) merge and transition from one to an- 
other. The connection between T = degenerate and 
Coulomb dominated objec ts has been known from the 
iZapolskv fc SalpeteJ l)1969j) M-R relations. The connec- 
tion between polytropes and degenerate objects, neglect- 
ing Coulomb physics, is seen in the n — 3/2 polytropes. 
Our models make the final connection between the three 
classes of objects, filling in the gap occupied by T 7^ 
objects in which Coulomb physics cannot be neglected. 

As discussed in §??, a ubiquitous feature of stellar M- 
R relations at sufficiently high Tc is the existence of 
a minimum mass, Mmin, below which equilibrium so- 
lutions do not exist. The cause of this is the transi- 
tion from ideal gas to degenerate electrons providing the 
pressure support and we showed in i )2.2l that the well 



known max imum in n = 3/2 polytropes fo r a fixed M 
llCox fc Gu iU 1968 : iRappaport fc JosjlTosl: iSteyensonl 
119911: iBurr ows fc Lieberti 119931: lUshomirskv et al.lll99'^ 



is the same physics. The value of Mmin depends on both 
Tc (in fact, Mmin = at sufficiently low Tc) and the 
strength of Coulomb physics. In general, the lower Tc 
and the stronger the Coulomb interactions, the smaller 
Afmin- The existence of an Mmin may have a profound 
impact on the evolution of a donor undergoing mass loss. 
For our He WDs, M„ii^ > for m odels with Tc > 10^ 
K. As discussed in iBildst'enI l)2002f) . the existence of an 
AiTmin > leads to the possibility of disrupting the donor 
through mass loss. This could be accomplished through 
stable mass loss down to M2 — Mmin or via a mass trans- 
fer instability. The latter will occur if the expansion of 
the donor under mass loss exceeds that of the RL. The 
entropy input needed to cause this instability and the 
fate of the donor are the subject of future work. 

We have applied our model to the accreting ultracom- 
pact MSPs. In XTE J 175 1-305, the donor must be hot re- 
gardless of its composition; in XTE J0929-314 and XTE 
J1807-294, fully degenerate donors are possible, depend- 
ing on the composition. From orbital inclination con- 
straints, the probabihty that XTE J1807-294 can ac- 
commodate a He donor is 15% while for XTE J0929- 
304, this probability is « 35%, providing support to the 
notion that some of these donors are likely C/0 WDs 
llSchulz et alJl2001t IjHett. Psaltis. fc Chakrabartvll200H 
iHomer eT^hl 200^ . The evolution of each system will 
differ depending on the orbital inclination. In particu- 
lar, how far the system can evolve in Porb in a specified 
time and the expected M-Porb relation depends on sini 
through both the mass, core temperature, and compo- 
sition of the donor. In general, the M-Porb relation is 
additionally parameterized by the donor Tc and composi- 
tion. We find that the number distribution of systems as 
a function of Porb, -^(Porb), is determined by the donor's 
riR. The distribution for systems with C/0 donors thus 
varies significantly from those with He donors. In the 
case of adiabatic evolution, the relative number of sys- 
tems at 40 min to that at 10 min is 160 for C/0 
donors and f» 60 for He donors. In addition to the ac- 
creting MSP systems, our models are applicable to the 
AM CVn binary systems, which ar e believed to b e dou- 
ble WD binaries with a He donor The 
application of our models to these systems is the subject 
of current work. 

The constant entropy models we have calculated give 
a lower limit on R for a given M, Tc, and composition. 
The actual thermal profile of a donor will depend on its 
prior evolution: how entropy is deposited into the star, 
entropy losses, and how quickly heat transport occurs as 
compared to mass loss. To determine the entropy profile 
of a donor in a specific system requires consideration of 
the coupled evolution of the binary and the donor. While 
our models do not address this uncertainty, they do pro- 
vide limiting M-R relations based on the total entropy 
content of the model in a consistent and systematic treat- 
ment without consideration of past evolution. As such 
they will be useful in consideration of binary systems 
with low-mass WD companions where time evolution of 
the donor's structure coupled to that of the binary it- 
self is not computationally feasible or necessary. Tracks 
of our M-R relations and rzad as a function of pc and 
Tc for He, C, and O donors are available along with the 
electronic version of the paper. 
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